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摘 E: KER (Hibiscus mutabilis) 栽培 历史 悠久 ,是 原 产 中 国 的 古老 园林 树种 和 药 用 植物 。 为 探讨 木 英 
车 品 种 及 近 缘 种 的 进化 特征 ， 厘 清 木 芙蓉 品种 间 及 其 与 近 缘 种 间 的 亲缘 关系 ， 同 时 也 为 探究 木 芙 蓉 叶 绿 体 
基因 组 (chloroplast DNA, cp DNA) 的 遗传 方式 , 该 文选 择 了 一 个 杂交 组 合 中 的 3 ARERR mP (“ 单 
WA . “金秋 人 掺 : 、“ 牡 丹 粉 ” ), 用 高 通 量 测序 平台 lumina NovaSeq 对 其 叶绿体 基因 组 进行 首次 测序 。 

经 组 装 注释 后 得 到 3 条 完整 的 cp DNA 序列 ， 结 合 本 团队 已 经 完成 的 近 缘 种 台湾 芙蓉 (所 taiwanensis) 和 来 
自 基因 库 的 木 樟 、 朱 栓 的 cp DNA, iZOOS AERIS 4 种 及 木 芮 敬 种 下 的 3 个 品种 进行 了 cp DNA 组 成 和 结 
构 特 征 的 比较 分 析 ， 并 完成 了 其 系统 发 育 树 重 建 。 结 果 表 明 : (0) HWA + “KS” 、“ 牡 丹 粉 ”的 
叶绿体 基因 组 序列 长 度 分 别 为 160 880. 160 879. 160 920 bp， 基 因数 目 均 为 130 个 ， 其 中 85 个 蛋白 编码 
基因 、8 个 核糖 体 RNA 和 37 个 转运 RNA。(2) 比 较 分 析 表 明 ， 木 黄蓉 的 种 下 3 个 品种 及 其 近 缘 种 台湾 芙蓉 
在 叶绿体 基因 组 上 高 度 保 守 , 反 向 重复 区 (IR) 均 为 26 300 bp; 木 桂 和 朱林 在 IR 区 发 生 了 收缩 ,分 别 为 25 745 
和 25 598 bp。(3) 系 统 发 育 分 析 表 明 ， 种 下 3 个 品种 聚 成 一 个 单 系 支 ， 再 与 台湾 芙蓉 聚 成 一 个 高 支持 率 分 
支 ， 表 明 木 芙蓉 和 台湾 芙蓉 的 杀 缘 关系 最 近 ; 相 较 于 木村 和 朱 梭 ， 木 芙蓉 、 人 台湾 甘 共 两 种 与 海滨 木 模 、 黄 
柳 、 大 麻 检 在 亲缘 关系 上 更 近 。(4) 木 美 蓉 3 个 品种 之 间 能 通过 叶绿体 基因 组 序列 区 分 开 ， 在 大 /小 单 拷贝 
区 (LSC/SSC) 长 度 上 ， HIWA n “金秋 颂 ”、“ 特 丹 粉 ”分 别 为 89 355 bp/18 925 bp. 89 353 bp/18 926 
bp、89 400 bp/18 920 bp， 且 从 重复 序列 和 核 苷 酸 多 样 性 分 析 中 开发 出 了 候选 分 子 标记 和 DNA 条 形 码 ， 可 
以 作为 品种 鉴定 的 分 子 条 码 。(5) 木 芙蓉 品种 单 兴 白 与 金秋 颁 叶 绿 体 基因 组 差异 最 小 ， 亲 缘 关 系 最 近 ， 根 据 
两 者 母 本 与 子 代 的 关系 ， 证 明了 木 芙 荧 叶 绿 体 基因 组 的 母系 遗传 特征 。 该 文 有 助 于 我 们 更 好 地 了 解 3 个 木 
芙 敬 品 种 及 台湾 芙蓉 叶绿体 基因 组 的 进化 特征 和 物种 间 的 系统 发 育 关 系 ， 并 为 木 芙蓉 品种 的 准确 鉴定 和 优 
良品 种 选 育 提供 了 叶绿体 基因 组 方面 的 基础 资料 。 
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Comparative analysis of the chloroplast genomes of three cultivars of 


Hibiscus mutabilis and its related species 


LI Zhenbing!, REN ting!, DENG Jiaojiao', CHEN Junpei', ZHOU Songdong", ZENG Xinmei’, 
MA Jiao’, LI Fangwen? 

(1. College of Life Sciences, Sichuan University, Key Laboratory of Bio-Resource and Eco-Environment of 
Ministry of Education, Chengdu 610065, China; 2. Chengdu Institute of Landscape Architecture, Chengdu 610083, 
China) 

Abstract: Hibiscus mutabilis is native to China with a long cultivation history, and is an ancient garden tree 
species and medicinal plant. In this study, we selected three cultivars of H. mutaibilis in a hybrid combination (H. 
mutabilis cv. Danbanbai, H. mutabilis cv. Jingiusong, H. mutabilis cv. Mudanfen) to investigate evolutionary 
characteristics between the cultivars of Hibiscus and its related species, and clarify the phylogenetic relationship 
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between the cultivars of Hibiscus and with its related species, as well as explore the genetic model of cp DNA of 
H. mutaibilis at the same time. We first sequenced the three cultivars of H. mutaibilis using Illumina NovaSeq. 
After assembly and annotation, three complete chloroplast genome sequences were obtained. The cp DNA of the 
related species H. taiwanensis from our group, and H. syricus and H. rosa-sinensis from the gene bank. Then we 
carried out comparative analysis on composition and structure of cp DNA of 4 species of Hibiscus and 3 cultivars 
of H. mutabilis, and completed its phylogenetic tree reconstruction. The results were as follows: (1) Total size of 
chloroplast genomes of H. mutabilis cv. Danbanbai, H. mutabilis cv. Jinqiusong, H. mutabilis cv. Mudanfen was 
160 880, 160 879, 160 920 bp, respectively, and the total gene number was 130, including 85 protein-coding 
genes, eight ribosomal RNAs, and 37 transfer RNAs. (2) The comparative analyses showed that the cp DNA of 
three cultivars of H. mutabilis and the related species H. taiwanensis were highly conserved, and the inverted 
repeat regions (IR) were all 26 300 bp; H. rosa-sinensis and H. syriacus shrank to at 25 745 and 25 598 bp, 
respectively. (3) The phylogenetic analysis revealed that the three cultivars were planted into a monophyletic 
branch, and then together with H. taiwanensis into a high support branch, indicating that H. mutabilis and H. 
taiwanensis had the closest relationship; Compared with H. syriacus and H. rosa-sinensis, H. mutabilis and H. 
taiwanensis were more closely related to H. hamabo, H. tiliaceum and H. canabinus. (4) Three cultivars of H. 
mutabilis could be distinguished by chloroplast genome sequence, the length of LSC/SSC of H. mutabilis cv. 
Danbanbai, H. mutabilis cv. Jinqiusong, H. mutabilis cv. Mudanfen were 89 355 bp/18 925 bp, 89 353 bp/18 926 
bp, 89 400 bp/18 920 bp, respectively. And candidate molecular markers and DNA barcodes had been developed 
from repeat sequence and nucleotide diversity analyses, which could be used as a tool for cultivars identification. 
(5) The chloroplast genomes of H. mutabilis cv. Danbanbai and H. mutabilis cv. Jinqiusong showed a minimum 
difference and had the closest phylogenetic relationship. According to the relationship between their female and 
offspring, the maternal genetic characteristics of the chloroplast genome of Hibiscus were proved. This study will 
help us to understand the evolutionary characteristics and phylogenetic relationship of chloroplast genomes of 
three cultivars of H. mutabilis and H. taiwanensis, and provide basic data on chloroplast genome for accurate 
identification of the cultivars of H. mutabilis and breeding of excellent cultivars. 

Keywords: Hibiscus mutabilis, Hibiscus taiwanensis, chloroplast genome, molecular marker, phylogenetic 
relationship 


AX (Hibiscus mutabilis) PERTE, ims RBS. BTE, CERE RAKIA 
布 。 其 花 大 色 艳 且 花 期 较 长 WEE, 2019). KS CRITIC. FA KEREJ, EIT p 
中 具有 重要 的 应 用 价值 ( 郑 觅 等 , 2012)。 成 都 市 植物 园 长 期 以 来 致力 于 市 花木 英 芍 的 研究 及 应 用 ， 并 培育 
了 许多 抗 病 抗 虫 品种 和 不 同 花 期 的 品种 ( 王 莹 ，2017)。 但 长 期 以 来 的 人 工 选 育 和 自然 杂交 导致 木 芙 蓉 品 种 
间 遗 传 关系 错综复杂 、 亲缘 关系 不 明 、 分 类 和 进化 关系 模糊 。 对 木 蔷 获 进行 品种 分 类 以 及 探讨 其 亲缘 关系 ， 
对 于 木 芙蓉 品种 间 杂 交 以 及 新 品种 的 选 育 具 有 重要 意义 OKRE, 2021). GS (Hibiscus taiwanensis) 
别名 狗头 黄蓉 或 山 芙蓉 ， 亦 属于 木 樟 属 ， 原 产 于 台湾 阿里 山 ， 是 当地 特有 的 原生 种 (Lim, 2014). 4h EE 
物 志 记载 ， 台 湾 芙 欧 非 常 类 似 于 木 芮 其， 区 别 仅 为 前 者 全 株 具 糙 硬 毛 和 粗糙 的 毛 被 ， 而 后 者 全 株 为 星 状 绒 
E, 因此 , 有 学 者 提出 两 者 可 能 不 是 明确 的 两 个 种 QI ER, 1984)。 此 外 , 木 芙蓉 与 台湾 芙 芍 很 容易 杂交 ， 
坐果 率 为 57.89%, 这 与 木 芙 蓉 品 种 之 间 杂 交 的 坐果 率 62.5% 相 差 无 几 , TU AN D 53 AER A AC AI TRES 7g 8.33%， 
与 木 樟 杂 交 成 功率 更 低 ， 这 在 某 种 程度 上 也 反映 了 亲缘 关系 的 远近 。 遗 传 信 息 不 足 通常 使 我 们 无 法 充分 了 
解 栽培 植物 和 近 缘 种 ， 明 确 栽 培 种 和 野生 近 缘 种 之 间 的 遗传 变异 对 将 近 缘 种 的 有 利 性 状 引入 栽培 品种 至 关 
重要 (Amar et al., 2019)。 因 此 ， 我 们 需要 开发 更 多 的 遗传 信息 ， 以 便 能 够 比较 栽培 种 的 遗传 差异 ， 并 快速 
和 准确 地 鉴定 它们 ， 从 而 更 有 效 地 利用 这 些 物种 。 

叶绿体 基因 组 (cp DNA) 大 小 为 75 一 250 kb, 大 多 数 植物 的 cp DNA 由 大 约 120 个 基因 组 成 , 包括 转运 
RNA (RNA)、 核 糖 体 RNA (rRNA) 和 香 白 质 编 码 基 因 (protein-coding gene, PCGs)(Daniell et al., 2016). cp 
DNA 因 结 构 简 单 、 高 度 保守 、 拷 贝 数 多 而 分 子 量 少 的 特点 ， 已 被 应 用 于 分 子 标 记 开 发 和 系统 发 育 等 研究 
中 ( 杨 俏 俏 等 ，2019)。 越 来 越 多 的 研究 表明 叶绿体 基因 工程 在 植物 遗传 改良 方面 具有 很 大 优势 ， 叶 绿 体 也 
己 经 成 为 植物 遗传 转化 的 新 工具 ( 母 连 胜 等 ，2017)。 前 人 关于 木 芙 蓉 的 研究 大 多 集中 在 木 美 荣 栽 培育 种 ( 王 
莹 ，2017; 石 小 庆 等 ，2021) 和 化 学 成 分 上 ( 蔡 露 等 ，2021; 王 艺 等 ，2021)， 鲜 少 有 关于 木 芙 蓉 品 种 间 以 及 
与 其 近 缘 种 间 杀 缘 关 系 的 研究 。 近 年 来 ， 有 学 者 在 木 芙蓉 和 人 台湾 芙蓉 的 叶绿体 基因 组 分 析 和 杀 缘 关系 研究 
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上 取得 了 一 定 的 进展 ，Abdullah 等 (2021) 第 一 次 对 包括 木 美菱 在 内 的 三 个 物种 的 cp DNA 进行 了 测序 和 比 
较 分 析 ， 但 其 研究 对 象 为 整个 锦 蒙 科 的 3 个 不 同属 植物 ， 范 围 偏 大 。 本 团队 的 许 新 锐 等 (2020) 对 台湾 其 
蕉 的 叶绿体 基因 组 进行 了 报道 ， 张 璐 等 (2021) 采 用 扫描 电镜 观察 了 一 些 木 芙蓉 品种 的 花粉 显 微 结 构 ， 并 探 
讨 了 其 分 类 学 意义 。 花 粉 稳定 性 极 强 ， 在 一 定 程 度 上 能 反映 植物 演化 和 亲缘 关系 ， 并 对 物种 鉴定 有 一 定 的 
分 类 价值 ( 彭 焕 文 等 ，2018)， 但 基因 组 层面 的 分 子 数 据 具 有 海量 的 遗传 信息 ， 通 常 认为 与 基因 组 层面 的 研 
究 相 结合 有 助 于 更 精确 和 深入 地 进行 亲缘 关系 和 物种 鉴定 的 研究 。 而 至 今 为 止 公布 的 木 芙 蓉 品 种 的 遗传 信 
息 很 少 ， 且 未 有 对 木 芮 蕉 品种 间 及 其 与 近 缘 种 台湾 鞭 鞭 间 的 cp DNA 相关 研究 。 这 就 极 大 阻碍 了 木 芙 蓉 品 
种 的 系统 发 育 和 优良 性 状 改良 研究 。 由 于 cp DNA 并 非 都 是 母系 遗传 ， 只 是 在 大 多 数 被 子 植物 中 为 母系 遗 
传 ， 在 裸子 植物 中 则 主要 为 父系 遗传 ， 此 外 ， 少 部 分 被 子 植物 的 cp DNA 存在 父系 遗传 和 双亲 遗传 的 现象 
(Neale & Sederoff, 1989), FAX Y TERZA BJ cp DNA 遗传 方式 ， 我 们 选取 一 个 杂交 组 合 中 的 3 ARR 
蕉 品种 ， 首 次 对 木 芮 蓉 品 种 及 其 近 缘 种 台湾 芙蓉 叶绿体 基因 组 进行 比较 分 析 和 系统 发 育 分 析 。 

本 研究 拟 探讨 以 下 科学 问题 : (1) 木 芙蓉 3 个 品种 及 其 近 缘 种 台湾 芙蓉 叶绿体 基因 组 有 什么 进化 特征 ? 
(2) 木 芙蓉 与 其 近 缘 种 间 有 怎样 的 亲缘 关系 ? (3) 在 叶绿体 基因 组 的 组 成 和 结构 方面 , 能 否 开发 出 品种 鉴定 的 
分 子 标 记 或 DNA 条 形 码 ? (4) 探究 木 芙 蓉 cp DNA 的 遗传 方式 。 本 研究 将 为 木 英 蕉 的 品种 鉴定 、 进 化 研究 
和 植物 遗传 改良 、 优 良品 种 选 育 提 供 重 要 的 第 一 手 遗 传 资 料 。 


1 材料 与 方法 


1. 1 木 芙 蓉 品 种 材料 来 源 

3 个 木 英 蓉 品 种 采 自 成 都 市 植物 园 (104*8'11" E. 30°45'52" N), DIA HIRA’ (H. mutabilis cv. 
Danbanbai)、“ 金 秋 颁 *(H. mutaibilis cv. Jinqiusong). EPI (H. mutabilis cv. Mudanfen)。 这 三 个 品种 为 杂交 
组 合 ， 其 中 ' 单 六 白 * 是 母 本 , “牡丹 粉 ' 是 父 本 ，' 金 秋 颂 ' 为 子 一 代 ， 这 样 有 利于 对 木 美 著 叶 绿 体 的 遗传 方式 
进行 探讨 。 从 生态 特性 来 看 ，' 单 淮 白 * 和 “牡丹 粉 ;为 早 花 品种 ， 花 期 为 6 一 9 月 ，'“ 金 秋 颁 为 中 花 品种 ， 花 
期 为 9 一 10 月 。 从 形态 上 来 看 ，'“ 金 秋 颁 "为 重 兴 花 ， 雄 蕊 退化 ， 有 子 房 但 不 结实 ， 另外 两 种 结实 。' 单 办 白 ” 
为 单 淮 花 ， 花 白色 ，' 特 丹 粉 ' 和 “金秋 颁 ' 分 别 为 粉色 和 红色 。' 牡 丹 粉 ?的 花 型 似 牡丹 ， 且 平均 花 径 比 前 两 者 
大 。 近 缘 种 台湾 芙蓉 的 花期 始 于 10 月 下 旬 ， 结 实 ， 花 为 单 闪 ， 白 中 带 粉 。 花 的 形态 见 图 版 I 所 示 。 
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a. HIA’ b “牡丹 粉 ”; c “金秋 颂 ”; d. SEXE. 
a. H. mutabilis cv. Danbanbai; b. H. mutabilis cv. Mudanfen; c. H. mutabilis cv. Jinqiusong; d. H. taiwanensis. 
图 版 1 3 个 木 芙 蓉 品 种 和 台湾 甘 蓉 花 的 形态 特征 
Plate I Morphological characters of flowers of three cultivars of Hibisucs mutabilis and H. taiwanensis 


1. 2 基因 组 DNA 的 提取 和 测序 


采集 3 ARRA nn 


种 新 鲜 健康 的 叶片 ， 用 改良 CTAB 法 从 叶片 组 织 中 提取 总 基因 组 DNA。 使 用 
脂 糖 凝 胶 电泳 和 荧光 染料 (QuantiT PicoGreen dsDNA Assay Kit) 评 估 DNA 产物 的 质量 和 浓度 。 
Illumina TruSeq Nano DNA LT 文库 制备 实验 流程 构建 插入 片段 大 小 (insert size) 为 400 bp 的 DNA 文库 。 


1% 琼 
采用 


EN 
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用 Ilumina NovaSeq 平台 对 DNA 文库 进行 二 代 测 序 (paired-end, 2x150 bp). DNA 提取 及 测序 工作 均 在 南京 
派 森 诺 基 因 科技 有 限 公司 完 成 。 
1.3 叶绿体 基因 组 的 组 装 和 注释 


每 个 物种 得 到 至 少 5G 的 原始 数据 (raw data), 经 过 数 
al, 2017) 软 件 (k-mer = 39) 进行 组 装 拼 接 ， YETE GBA 
利用 在 线程 序 GeSeq(Tillich et al., 2017) 注 释 叶 绿 体 基 
行 手动 校正 ， 最 后 将 带 有 注释 信息 的 3 条 序列 上 传 至 NCBI Z 


因 组 序列 ， 并 用 


BÆ, “Ñ 


La 


居 过 滤 后 ， 首 先 使 用 NOVOPlasty (Dierckxsens et 
的 rbcL 基因 作为 种 子 序列 (seed sequence)。 然 后 
Geneious v9.0.2(Kearse et al., 2012) 进 


“金秋 颂 ) “牡丹 粉 " 的 序列 


号 分 别 为 MZ846191、MZ846192、MZ855502。 借 助 在 线程 序 Oranellar Genome DRAW(Greiner et al., 2019) 


绘制 叶绿体 基因 组 的 物理 图 谱 。 
本 文中 用 到 的 其 他 叶绿体 基因 组 序列 均 下 载 于 NCBI, 
1.4 叶绿体 基因 组 比较 分 析 


台湾 芙蓉 的 序列 号 为 MK937807,， 亦 为 本 研究 项 目的 成 果 (CXu et al., 2019). 


朱 检 基因 组 


REPuter(Kurtz et 


lim 


置 为 3， 重 复 序 列 的 最 小 长 度 限 制 为 30)。 借 助 IRscope (Amiryousefi et al., 2018)22 EIN ER 3 个 品种 、 台 
组 IR 边界 的 扩张 和 收缩 。 使 用 mVISTA(Mayor et al.，2000) 中 的 
木林 、 朱 樟 的 叶绿体 基因 组 序列 进行 全 局 比 对 ， 以 
E. 并 用 DNAsp6(Librado & Rozas，2009) 软 件 计 算 序 列 之 间 的 核 背 酸 多 态 性 。 使 用 
Geneious v9.02(Kearse et al., 2012) 分 别提 取 每 个 物种 的 和 蛋白 编码 区 域 (CDS)， 存 在 多 拷贝 的 基因 仅 保 留 一 
然后 用 Codon W(Sharp & Li, 1987; £ 


zb Ih 


BRER, 


SE 


usage, RSCU). 
1.5 系统 发 育 分 析 


后 用 concatenate 功能 将 序列 是 


X Bë JI K FE i: 2 s dis [X| 
Shuffle-LANGAN 模式 对 木 芙蓉 3 个 品种 、 
评估 序列 之 间 的 相似 怡 


使 用 phylosuite(Zhang et al., 2020) 
的 基因 后 ， 用 MAFFT 对 序列 进行 比 对 ， 而 后 


=n 
BER. f 


使 用 MISA(MicroSatellite identification tool)(Beier et al., 2017) 15 
中 的 微 


ANE 


r1 SA. 


提取 17 MHF 


] Modeltest 软件 检 洲 


ita 


| 信息 见 表 1。 


这 些 CDSs 的 同 义 密码 子 的 相对 使 ) 


IKE 3 Amih GEER, 
gH (SSRs)。 参 数 设 置 : 单 核 苷 酸 (mono)、 双 核 
Hë (tetras MIZER (penta), NIZE (hexa) SSR 的 最 小 习 


al., 2001) 来 分 析 木 芙蓉 3 个 品种 和 人 台 


HE? (di), ZER (tri). MRE 
E UM 10. 5. 4. 3. 3. 3. AIH 
E KAIRE AE rH I 8 E 2S7 (Hamming 距离 设 


zi 十 ht 


Aë. 


条 , 


叶绿体 基因 组 的 CDSs， 删 去 重复 基因 和 非 共 有 
] MACSE 优化 序列 ， 接 着 用 Gblock 修剪 比 对 好 的 序列 ， 最 
上 最 佳 核 苷 酸 蔡 代 模 型 为 JC+IHG。 使 用 


JÆ (related synonymous codon 


MrBayes 


软件 进行 贝 叶 斯 (bayesian inference, BD 分 析 。 基 于 Markov chain Monte Carlo MCMC) 算 法 , 运行 一 千 万 代 ， 
每 1000 代 取 一 棵 树 ， 前 20% 的 树 当 作 老 化 样本 丢弃 ， 
然 法 (maximum likelihood, ML) 分 析 , 重复 次 数 为 1 000。 最 后 将 建树 结果 用 在 线程 序 iTOL(Ivica & Peer, 2021) 


剩余 的 树 ) 


j 于 构建 一 致 树 。 使 用 


IQtree 进行 最 大 似 


H. tiliaceum 


Malva verticillata 


进行 美化 。 
表 1 本 研究 所 选用 的 物种 名 称 和 GenBank 序列 编号 

Table 1 Names and GenBank sequence numbers of the plant species selected in this study 
物种 序列 编号 物种 序列 编号 
Species Accession number Species Accession number 
AE KR259989 桐 棉 NC048518 
Hibiscus syriacus Thespesia populnea 
ATE NC042239 苘 麻 NC053702 
H. rosa-sinensis Abutilon theophrasti 
AREE NC045873 DE NC053839 
H. canabinus Alcea rosea 
EE MT644160 Ep MT106775 
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海滨 木 樟 NC030195 ¿HEH NC051877 
H. hamabo Sida szechuensis 
海岛 棉 NC008641 美丽 异 木棉 MK820674 
Gossypium barbadense Ceiba speciosa 
G. thurberi GU907100 

2 结果 与 分 析 


2.1 叶绿体 基因 组 基本 特征 
如 图 1 所 示 ，“ 单 瓣 白 ” 和 “金秋 颂 ” 的 长 度 差异 较 小 ， 分 别 为 160 880 和 160 879 bp, 
碱 基 ， 而 “牡丹 粉 ”的 长 度 与 前 两 种 则 有 稍 大 的 差异 ， 为 160 920 bp。 三 者 的 cp DNA 呈现 


仅 相 差 一 个 
典型 的 四 分 


体 环 状 结构 ， 均 由 一 对 反 辣 重复 区 (IRa 和 IRb， 长 度 均 为 26 300 bp) 以 及 一 个 大 单 拷贝 区 (LSC， 长 度 分 
别 为 89 355. 89 353、89 400 bp) 和 一 个 小 单 拷贝 区 (SSC, 长 度 分 别 为 18 925 bp. 18 926 bp. 18 920 bp) 组 


Ko AR 2 可 以 看 出 台湾 芙 苏 的 LSC 区 长 度 亦 明显 长 于 3 个 木 黄蓉 品种 ， 而 在 品种 内 部 ， 
“金秋 颂 ” 的 LSC 区 与 SSC 区 长 度 仅 相差 1 一 2 bp， 但 与 “牡丹 粉 ” 相 差 较 大 。 
木 黄 蕉 的 叶绿体 基因 组 一 共有 130 个 基因 ,包含 85$ 个 蛋白 编码 基因 (PCGs)、37 个 tRNA 和 
大 多 数 基因 以 单 拷贝 形式 出 现 于 LSC 区 或 SSC 区 中 ,其 中 SSC 区 中 有 13 个 基因 ,包括 12 个 P 


“CORE 5 


8 个 IRNA。 
CGs (ndhF. 


rpl32、 ccsA、, ndhD. psaC, ndhE, ndhG. ndhl, ndhA, ndhH. rps15、 ycfl) 和 1 个 tRNA (trnL-UAG), LSC 
区 中 有 85 个 基因 , 包括 63 个 PCGs 和 22 个 tRNA。 NA 17 个 基因 在 IR 区 重复 , 包括 6 个 PCGs (rpl2 , 
rpl23 ,ycPP .ndhB\rps7 .rps12).7 个 tRNA (trnI-CAU trnL-CAA \trnV-GAC , trnI-GAU trnA-UGC ,trnR-ACG. 
trmnN-GUU) 和 4 个 IRNA (rm16, rm23, rrn4.5, rrn5). ycfl 基因 跨越 SSC 区 和 IR X, rps12 基因 第 一 个 外 
显 子 位 于 LSC 区 , 其 他 两 个 外 显 子 位 于 IRs 区 。17 个 基因 有 一 个 内 含 子 , ycf3 和 clpP 基因 有 两 个 内 含 子 。 


所 有 基因 按照 基本 功能 划分 为 表达 相关 基因 、 光 合作 用 相关 基因 和 其 他 基因 三 类 ， 每 类 基因 的 名 称 和 数量 


展示 在 表 3 中 。 


HRH’ 
H mutabilis cv. Danbanbai 
160 880 bp 


“金秋 颁 ” 
H. mutabilis cv. Jinqiusong 
160 879 bp 


“牡丹 粉 ” 
H. mutabilis cv. Mudanfen 
160 920 bp 


mi Ho Hs Bh Ds M6 M7 Els Mo Mio Bs Mi: 国 13 口 14 


| 


四 


1. 光合 系统 I; 2. 光合 系统 II; 3. 胞 色素 复合 物 ; 4. ATP. 合成 酶 ，5. NADH MAR 6. 核 酮 糖 二 磷酸 羧 化 酶 大 亚 基 ; 7. RNA 


聚合 酶 ，8. 核糖 体 蛋 白 小 亚 基 (SSU); 9. 核糖 体 蛋 白 大 亚 基 (LSU); 10. 转运 RNAs; 11. 核糖 体 RNAs; 12. 


clpP. matK 2= 


D]; 13. 其 他 基因 ; 14. 假设 的 叶绿体 阅读 框 Qyc/)。 
1. Photosystem I; 2. Photosystem II; 3. Cytochrome b/f complex; 4. ATP synthase; 5. NADH dehydrogenase; 6. RubisCO large 
subunit; 7. RNA polymerase; 8. Ribosomal proteins (SSU); 9. Ribosomal proteins (LSU); 10. Transfer RNAs; 11. Ribosomal RNAs; 
12. clpP, matK; 13. Other genes; 14. Hypothetical chloroplast reading flames (ycf). 

图 1 木 芙 蓉 叶 绿 体 基 因 组 物理 图 谱 
Fig. 1 Physical map of Hibiscus mutabilis chloroplast 


K 2 EE ARE BERRA 3 个 品种 ) 叶绿体 基因 组 基本 特征 比较 


Table 2 Comparison of chloroplast genomes among four species of Hibiscus (including three cultivars of H. 


mutabilis) 
KER KÆ Size (bp)/ GC 含量 GC content (%) 基因 数 
Hibiscus Number of gene 
全 基因 组 ”小 单 拷贝 区 大 单 拷贝 区 反 向 重复 区 总 和 蛋白 编码 基因 ”核糖 体 转运 
Genome SSC LSC IR Total PCGs RNA/rRNA RNA/RRNA 
HA’ 160 880/36.9 18 925/31.5 89 355/34.7 26 300/42.6 130 85 8 37 
H. mutabilis cv. Danbanbai 
“AKAN 160 879/36.9 18 926/31.5 89 353/34.7 26300/42.6 130 85 8 37 
H. mutabilis cv. Jinqiusong 
“ 特 丹 粉 ” 160 920/36.9 18 920/31.5 89 400/34.7 26 300/42.6 130 85 8 37 
H. mutabilis cv. Mudanfen 
台湾 黄蓉 161056/36.9 18918/31.5 89 538/34.7 26300/42.6 130 85 8 37 
H. taiwanensis 
AE 160 951/37.0 20246/31.3 89 509/34.9 25 598/42.9 130 85 8 37 
H. rosa-sinensis 
AGE 161022/36.8 19831/31.1 89 701/34.7 25 745/42.8 130 85 8 37 


H. syriacus 


R 3 木 芙 蓉 叶 绿 体 基 因 组 的 基因 信息 


Table 3 List of genes in the chloroplast genome Hibiscus mutabilis 


基因 类 型 数量 
A SR 基因 分 组 Group of gene 基因 名 称 Name of gene 
Gene category Number 
rpl2a*, 23a, 32, 22, 16*, 14, 36, 
核糖 体 大 亚 基 基 因 Large subunit of ribosome gene B Iaa 11 
rps]a, 15, 19, 3, 8, 11, 12a*, 18, 
核糖 体 小 亚 基 基因 Small subunit of ribosome gene 14 
4, 14, 2, 16* 
RNA 聚合 酶 亚 基 基 因 
rpoA, rpoB, rpoCl*, rpoC2 4 


DNA-dependent RNA polymerase gene 

核糖 体 RNA 基因 Ribosomal RNA gene rrnl6a, rrn23a, rrn4.5a, rrn5a 8 
trnl-CAUa, trnL-CAAa, trnV-GACa , 
trnA-UGCa*, trnR-ACGa, trnI-GAUa*, 
trnN-GUUa, trnL-UAG, trnP-UGG , 
trnW-CCA, trnM-CAU, trnV-UAC* , 

转运 RNA 基因 Transfer RNA gene trnF-GAA , trnL-UAA* ,. trnT-UGU , 2 

trn$-GGA , trnfM-CAU , trnG-GCC ， 
trn$-UGA , trnT-GGU , trnE-UUC , 
trnY-GUA , trnD-GUC , trnC-GCA ， 
trnR-UCU , trnG-UCC*, trn$-GCU, 
trnQ-UUG, trnK-UUU*, trnH-GUG 


表达 相关 基 医 


Self-replication 
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光合 系统 I 基因 Photosystem I gene psaC, J, L A, B 5 
sbH, N, T, B. E, F, L, J, Z, G, 
光合 系统 [I 基因 Photosystem II gene á 15 
D, M, L K, A 
NADH 氧化 还 原 酶 亚 基 基 因 ndhBa*, H, A*, L G, E, D, F, C, 3 
Subunits of NADH-dehydrogenase gene K, J 
光合 作用 - 
细胞 色素 复合 物 亚 基 基 因 

Photosynthesis I petD*, B*, G, L, A, N 6 

Subunits of cytochromeb/f complex gene 
ATP 合 酶 亚 基 基 因 Subunits of ATP synthase gene atpB, E, I, H, F*, A 6 
核 酮 糖 二 磷酸 羧 化 酶 大 亚 基 基因 25: à 

rbc 
Large subunit of rubisco gene 
翻译 起 始 因 子 Translation initiation factor infA 1 
ATP 依赖 性 蛋白 酶 亚 基 p 基因 

. clpP** 1 

ATP-dependent protease subunit p gene 
成 熟 酶 基因 Maturase gene matK 1 
OREHA] Envelope membrane protein gene cemA 1 
其 他 基因 乙酰 辅酶 A 羧 化 酶 亚 基 基 因 T i 

acc. 
Other gene Subunit of acetyl-CoA-carboxylase gene 
c 型 细胞 色素 合成 基因 

ccsA 1 

C-type cytochrome synthesis gene gene 

假设 的 叶绿体 阅读 框 

ycf2a, ycfl, ycf4, ycf3** 5 


Hypothetical chloroplast reading frames(ycf) 

总 和 Total 130 
*, 一 个 内 含 子 ;，**. 两 个 内 含 子 ，a. 重复 基因 。 
*. One intron; **. Two introns; a. Duplicated gene. 
2.2 IR 边界 分 析 

通过 比较 木 贡 芍 3 个 品种 、 台 湾 黄蓉、 朱 梯 、 木 树叶 绿 体 基 因 组 的 IR/LSC 和 IR/SSC 边界 区 域 的 基因 

分 布 状况 , 识别 IR 的 扩张 或 收缩 。 如 图 2 所 示 , IR/LSC 和 IR/SSC 边界 附近 分 布 的 基因 包括 rps19 , rpl2 , 
ycfl、ndhF、trnH，, KRR 3 个 品种 与 台湾 英 蓉 之 间 IR 边界 一 致 。6 条 序列 的 SSC/IRa 边界 均 为 ycfl ÆA, 
在 木 芙蓉 和 台湾 芙蓉 SSC 区 长 度 均 为 4026 bp, TEE AERE] SSC 区 长 度 分 别 为 5599 和 5 083 bp. 
ndhF 基 因 在 6 条 序列 中 均 位 于 SSC 区 内 , 其 中 木 芙蓉 和 台湾 芙蓉 的 adjpF 基 因 距 SSC/IRb 边界 均 为 32 bp, 
其 余 两 种 距 边 界 150 bp. 同样 的 , rpI2 基因 在 6 条 序列 中 均 位 于 Rb 区 内 , 其 中 , KOS 16115 2 š IJ rp 人 2 
基因 距 LSC/IRb 边界 均 为 103 bp, AK TEE. 67 bp， 木 检 距 边界 113 bp. 


JLB JSB JSA JLA 


an p=ll3 a 150 bp. o 1 
m ET 7 ETT LN 
H. syriacus Immun m, D ——3À 


1 
130p 
K = isoto TM 


ABl ml22 / / 33795 ed wth) 


EI 32bp 


PETS m2 - EET T jos A Im 77] 
H. taiwanensis CETTE uu m, PEL ES —À 


bp 


» " 
“牡丹 粉 ? mp2 # m? 4 626 bp | L089 bp . 2P, nl 


H. mutabilis cv. Mudanfen r i a a a 


Ni 
" 29. MUI TI á š 
“KD m22 4626 bp |I 089 bp i 2bp , nH 


H. mutabilis ev. Jinqjusong [ER 生生 于) 下 


29bp ~ 103 bp. = 32 bp " š 
ET E p22 N 4626 bp J1 089 bp 


2bp y nH 
H. mutabilis ev. Dambanbai [a p mpg M — a ERE ——3À 
ji ` 103 bp. 


图 2 GEPON. AE, ARERECKORTE 3 个 品种 叶绿体 基因 组 的 反 向 重复 区 (GR)、 大 单 拷贝 区 (LSC). ^ 
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单 拷贝 区 (SSC) 边界 对 比 
Fig. 2 Comparation of the junction between LSC, SSC and IR regions of chloroplast genome among Hibiscus 
taiwanensis, H. syriacus, H. rosa-sinensis and three cultivars of H. mutabilis 


2.3 叶绿体 微 卫 星 分 析 和 重复 序列 分 析 

微 卫 星 序列 (microsatellite DNA)， 也 叫 简单 序列 重复 (simple sequence repeats，SSR)， 是 1 到 6bp 的 
重复 序列 ， 广 泛 分 布 于 叶绿体 基因 组 。SSR 具有 高 度 多 态 性 和 特异 性 ， 是 研究 基因 流 、 种 群 遗 传 学 和 基因 
作 图 的 有 价值 标记 (Chen et al., 2015)。 本 文中 , 我 们 分 析 了 6 种 SSR ( 单 、 二 、 三 、 四 、 五 、 六 核 昔 酸 SSR) 
在 6 条 叶绿体 基因 组 中 的 分 布 。 如 图 3 所 示 ， 数 量 最 多 的 均 为 单 核 音 酸 重复 ， 占 比 从 66.67% 到 74.55% 不 
等 ， 均 比 其 他 所 有 重复 类 型 加 起 来 还 多 。 木 芙蓉 3 个 品种 的 SSR 总 数 为 95 ， 台 湾 黄 芍 与 其 非常 近似 ， 为 
96。 作 为 对 照 ， 朱 橙 和 木 样 分 别 为 63 和 110， 差 异 较 大 。 木 黄蓉 3 个 品种 和 台湾 芙 芍 的 基因 组 中 存在 六 种 
完全 微 卫 星 序列 ， 而 朱 橙 和 木 樟 中 不 存在 六 核 苷 酸 重 复 。 木 芮 区 3 个 品种 与 台湾 芙蓉 的 共同 区 别 仅仅 在 于 
前 者 均 只 有 1 个 五 核 苷 酸 重复 ， 而 后 者 存在 2 个 五 核 苷 酸 重 复 。 木 芙蓉 3 个 品种 中 ，“ 牡 丹 粉 ”明显 区 别 
于 另外 两 种 ，“ 牡 丹 粉 ”的 单 核 音 酸 重 复 比 另外 两 种 多 一 个 ， 而 六 核 昔 酸 重复 比 另 外 两 种 少 一 个 。 


AB 


H. syriacus 


AR 


H. rosa-sinensis 


amd 


H. taiwanensis 


“牡丹 粉 ” 


H. mutabilis cv. Mudanfen 


“金秋 颂 ” 


H. mutabilis cv. Jinqiusong 
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H. mutabilis cv. Danbanbai 
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图 3 RSRK, AGE. ZRERECROE 3 个 品种 叶绿体 基因 组 的 微 卫星 分 析 (SSR) 


Fig. 3 SSR analysis of chloroplast genome in Hibiscus taiwanensis, H. rosa-sinensis, H. syriacus and three 


cultivars of H. mutabilis 

之 前 的 研究 表明 重复 的 出 现 表 示 该 位 点 是 基因 组 重 构 的 主要 热点 ， 这 些 重复 对 于 开发 种 群 遗传 学 研究 
中 的 遗传 标记 具有 重要 价值 (Bull et al., 1999)。 重 复 序 列 有 四 种 类 型 ， 即 正 向 重 (forward)、 反 癌 (reverse). 
互补 (complement)、 回 文 (palindromic)。 如 图 4 rs, KEX 3 个 品种 中 均 检测 到 四 种 类 型 的 重复 ， 而 台 
湾 芙 敬 无 互补 重复 ; KER 3 个 品种 均 为 22 个 正 向 重复 和 3 个 反 向 重复 ， 而 台湾 美 蓉 的 正 向 和 反 向 重复 
分 别 为 26 和 1。 木 芙蓉 3 个 品种 间 的 四 种 重复 的 比例 十 分 类 似 ， 但 同时 互相 之 间 还 是 存在 着 差异 。 虽 然 三 
者 的 正 向 重复 均 为 22 个 ， 但 其 构成 略 有 差别 ，“ 牡 丹 粉 ” 重复 单元 长 度 为 30—39 和 40 一 49 的 正 向 重复 
数量 分 别 为 16 和 6， 而 其 他 两 种 都 分 别 是 17 和 5; “RE GR AKAN 的 回 文 重复 均 为 12 个 ,而 “ 特 
丹 粉 * 为 11 个 ; ，“ 单 注 白 ” 和 “牡丹 粉 ”的 互补 重复 均 为 2 个， 而 “金秋 颁 ” 为 1 个 。 
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A.4 种 重复 序列 类 型 的 数量 ，B. 不 同 长 度 的 重复 序列 数量 。 
A. Number of four repeat types; B. Number of repeat sequences by length. 
图 4 台湾 芙蓉 和 木 芙 其 3 个 品种 叶绿体 基因 组 的 重复 序列 分 析 


Fig. 4 Analysis of repeated sequences of chloroplast genome in Hibiscus taiwanensis and three cultivars of H. 


mutabilis 


2.4 核 苷 酸 多 态 性 分 析 

利用 mVISTA 在 线程 序 ， 比 对 了 木 芙 蓉 3 个 品种 和 台湾 芙蓉、 朱 检 、 木 樟 的 叶绿体 基因 组 。 其 中 ，“ 单 
办 和 白 ” 作 为 参考 序列 ， 朱 樟 和 木 樟 则 作为 对 照 。 从 图 5 中 可 以 明显 观察 到 朱 樟 、 木 樟 均 与 参考 序列 差异 较 
大 ， 而 “牡丹 粉 ”、““ 金 秋 颂 ”、 台 湾 芙蓉 则 与 参考 序列 差异 较 小 。5 个 cp DNA 中 高 度 分 化 的 区 域 主 要 
位 于 基因 间隔 区 ， 但 一 些 和 蛋白 编码 区 如 ycfl 基因 变异 度 也 较 大 ， 而 最 保守 的 4 个 rRNA EA, ERER 3 
个 品种 和 台湾 芙蓉 中 完全 一 致 。 

为 了 更 加 精确 地 反映 木 甘 攻 3 个 品种 之 间 以 及 与 台湾 芙 蕉 间 的 核 背 酸 多 态 性 ， 我 们 分 别 计算 了 “金秋 
颂 ”、 “牡丹 粉 ”、 台 湾 美 蓉 的 CDS 和 基因 间隔 区 相对 于 参考 序列 单 汶 白 的 核 苷 酸 多 样 性 值 (nucleotide 
diversity values, Pi)。 如 图 6 所 示 , CDS 序列 较为 保守 ,其 Pi 值 普遍 较 小 。 仅 就 CDS 序列 来 看 ， “金秋 颁 ” 
仅 在 ycfl 基因 处 出 现 多 态 性 ，“ 牡 丹 粉 ” 仅 在 ndhB 基因 处 出 现 ， 台 湾 英 芍 在 accD. atpA, clpP., ndhB, 
ndhD、matK、ycfl 均 出 现 ， 但 Pi 值 最 大 也 不 超过 0.001 7。 而 在 基因 间隔 区 中 ， 台 湾 芮 营 有 14 个 基因 间隔 
区 有 差异 ， 其 中 trnR-UCU~atpA、psbZ~trnG-GCC、infA~rps8、ndhE~ndhG 这 四 个 基因 间隔 区 为 台湾 芙蓉 


E 


K, IR 区 中 Pi 值 较 小 ， 均 小 于 0.020 9。“ 特 丹 粉 ”有 3 个 基因 间隔 区 有 多 态 性 ， 其 中 psbZ-trnG-GCC. 
ycfA-cemA 为 牡丹 粉 的 热点 变异 区 域 ， 且 这 2 个 热点 都 在 LSC K. TA, R 区 的 变异 低 于 LSC 区 和 SSC 
Ko “金秋 颁 ” 无 基因 间隔 区 差异 。 这 些 变异 度 较 高 的 区 域 ， 可 用 于 设计 特定 的 DNA 条 形 码 。 
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Y 轴 表 示 介 于 50% 到 100% 之 间 平 均 的 相似 性 百分比 。 
The Y axis represents the percent identity ranging from 50% to 100%. 
图 5 6 个 叶绿体 基因 组 的 序列 相似 度 图 (以 单 欠 白 为 参考 ) 
Fig. 5 Sequence identity plot of the six chloroplast genomes(using Hibiscus mutabilis cv. Danbanbai as a 
reference) 
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X axis represents the names of protein-coding genes and genes intergenic regions; Y-axis represents nucleotide diversity of each 


window. 


图 6 RERA “RD . “KANM 5 “RÉ AERE TEET 
Fig. 6 Comparative analysis of nucleotide variability among Hibiscus taiwanensis, H. mutabilis cv. Mudanfen, H. 
mutabilis cv. Jinqiusong and H. mutabilis cv. Danbanbai 


2.5 选择 压力 分 析 和 密码 子 偏好 性 分 析 

相对 同 义 密码 子 使 用 度 (RSCU) 用 于 评估 编码 序列 中 同 义 密码 子 的 使 用 情况 , 其 中 RSCU 越 大 表示 偏 
好 性 越 强 。 我 们 选取 了 木 芮 蓉 3 个 品种 、 台 湾 芙蓉 以 及 其 他 13 个 相关 物种 ， 测 定 了 每 个 物种 叶绿体 基因 
组 的 密码 子 的 数量 以 及 RSCU。 计 算 结 果 表 明基 因 组 中 亮 氨 酸 含量 最 高 ， 异 亮 氨 酸 和 甘氨酸 其 次 ， 编 码 半 
胱 氨 酸 的 密码 子 数量 最 少 ， 次 少 为 色 氨 酸 和 丝氨酸 。 如 图 7 所 示 ， 除 色 和 氨 酸 和 甲 硫 氨 酸 外 ， 所 有 氨基 酸 都 
使 用 两 个 或 多 个 同 义 密码 子 ， 精 氨 酸 和 亮 氨 酸 均 由 6 个 同 义 密码 子 表 达 。RSCU 值 较 大 的 密码 子 大 部 分 为 
A 或 UU 结尾 , 即 基因 上 的 碱 基 A/T, 这 与 前 人 的 研究 结论 一 致 ， 即 植物 中 AM 结尾 的 密码 子 使 用 偏好 (Xie 
et al., 2018)。 如 图 7 所 示 ， 锦 葵 科 植物 的 密码 子 偏好 表现 出 较 高 的 保守 性 ， 但 不 同属 之 间 还 是 存在 着 一 定 
的 差异 。 木 美菱 3 个 品种 和 台湾 芙蓉 紧 紧 地 聚 在 一 起 ， 其 余 物 种 也 大 致 是 按照 属 和 族 进 行 聚 类 。 
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红色 和 白色 分 别 表 示 较 高 和 较 低 的 RSCU 值 。 图 的 右边 是 物种 间 的 系统 发 育 关 系 。 
Red and white indicate higher and lower RSCU values, respectively. The right of the figure shows the phylogenetic relationship 


among species. 
E]7 15 个 物种 (@ 3 个 品种 ) 所 有 叶绿体 基因 组 蛋白 编码 区 (CDS) 同 义 密码 子 的 相对 使 用 度 (RSCU) 
Fig. 7 RSCU values of all protein-coding (CDS) genes for chloroplast genomes from 15 species (including three 


cultivars) 

2.6 系统 发 育 关 系 

我 们 选取 16 个 锅 葵 科 植 物 的 叶绿体 基因 组 来 探讨 木 芮 芍 与 其 近 缘 种 间 的 杂 缘 关系 ， 并 以 木棉 科 植 物 
美丽 异 木棉 (Ceiba speciosa) 作为 外 类 群 。 使 用 17 个 叶绿体 基因 组 中 共有 的 76 个 CDSs 进行 最 大 似 然 分 
Jr (ML) 和 贝 叶 斯 分 析 (BD。 结 果 显 示 ，ML 和 BI 分 析 的 拓扑 结构 完全 相同 ， 并 且 绝 大 部 分 的 进化 支 都 
具有 很 高 的 后 验 概率 (posterior probabilities) 和 上 自 展 值 (bootstrap values)。 如 图 8 Hrs, KRX 3 个 品种 
和 台湾 芙蓉 同属 于 木 橙 属 ， 并 且 聚 在 同一 枝 亲 缘 关 系 最 近 ， 自 展 值 和 后 验 概率 均 为 最 大 值 。 
KPE H. cannabinus 
WE H. tiliaceum 
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" * [— — — — Gossypium thurberi 


海岛 棉 Gossypium barbadense 
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分 支 上 的 数字 分 别 代表 自 展 值 和 后 验 概率 。* 代表 两 种 分 析 为 最 大 支持 度 。 

The numbers on the branch are bootstrap support values and posterior probability. * represents maximum support in all two analyses. 
图 8 基于 ML 和 BI 分 析 构 建 的 15 个 物种 〈 含 3 个 品种 ) 共有 的 76 个 蛋白 编码 序列 系统 发 育 树 

Fig. 8 Phylogenetic relationships of the 15 species (including 3 cultivars) inferred from ML and BI analyses based 


NS 


on 76 Shared protein-coding genes 
3 讨论 与 结论 


3.1 木 芮 蔬 及 近 缘 种 的 叶绿体 基因 组 比较 

叶绿体 基因 组 比较 分 析 表 明 ， 木 芙 芍 的 种 下 3 品种 及 其 近 缘 种 台湾 芙 芍 在 cp DNA 上 高 度 保守 ， 反 癌 
重复 区 (IR) 均 为 26 300 bp， 木 芙蓉 与 该 属 的 木 樟 和 朱 检 比较 ， 后 二 者 在 IR 区 发 生 了 收缩 分 别 为 25 745 
和 25 598 bp. 3 KJ uf cp DNA 上 的 大 /小 单 拷贝 区 在 长 度 上 存在 差异 ， 且 “ 单 锥 白 ”和 “金秋 颂 ” 
差异 更 小 ， 在 大 /小 单 拷贝 区 (LSC/SSC) KEE, "BRE. “金秋 颁 ”、“ 牡 丹 粉 ”分 别 为 89 355bp/18 
925bp、89 353bp/18 926bp、89 400bp/18 920bp。GC 含量 是 判断 物种 间 杀 缘 关 系 的 重要 指标 (Li et al., 2018), 
KER 3 个 品种 与 台湾 芙蓉 的 总 GC 含量 、 耻 区 和 SC 区 的 GC EEEE ZKEERUR ERU GC 含量 
则 有 明显 差异 。 此 外 ， 木 芙蓉 3 个 品种 与 台湾 芙蓉 的 IR 区 边界 分 布 的 ycfl. trnH. ndhF. rpl2. rps19 基 
因 均 未 发 生 扩张 或 收缩 ， 而 同属 的 木 植 和 朱 栓 相对 木 芙 蓉 和 台湾 芙 芍 发 生 了 明显 的 波动 。 不 同 物种 之 间 的 
IR 边界 是 不 同 的 , 其 中 IR 边界 的 波动 是 造成 这 种 差异 的 主要 原因 (Ravi et al., 2008)。 因 此 ，GC 含量 和 IR 
区 边界 的 情况 在 较 大 程度 上 说 明了 木 芙蓉 与 台湾 芙蓉 的 亲缘 关系 非常 之 近 ， 而 木 桦 和 朱 栓 的 亲缘 关系 则 相 
对 较 远 。 
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同 义 密 码 子 是 由 突变 产生 的 ， 进 化 压力 导致 这 些 同 义 密 码 子 的 使 用 频率 发 生变 化 (Hanson & Coller, 
2017)。 密 码 子 偏好 性 是 物种 对 其 碱 基 组 成 、tRNA 丰 度 和 环境 选择 压力 长 期 适应 的 结果 (Novoa & de 
Pouplana, 2012; Behura & Severson, 2012)。 此 外 ， 密 码 子 偏好 性 会 影响 翻译 的 起 始 、 延 伸 和 准确 性 、mRNA 
的 前 切 和 和 蛋白质 的 折 闭 ( 任 桂 萍 等 ，2019)。 因 此 ， 密 码 子 偏好 性 在 一 定 程度 上 也 能 反映 出 亲缘 关系 。 文 中 
密码 子 偏好 性 聚 类 结果 将 台湾 芙蓉 与 木 芙蓉 3 个 品种 聚 在 一 起 ，“ 金 秋 颁 ”与 “ 单 办 白 ”的 密码 子 偏好 性 
无 差异 ，“ 牡 丹 粉 ” 和 台湾 甘 鞭 相 比 前 者 编码 亮 氨 酸 密码 子 UUA Mfr MENSES, Hohe SAM TERR 
酸 和 丝氨酸 的 密码 子 偏 好 上 表现 出 差异 ， 所 以 从 密码 子 偏 好 性 上 我 们 也 能 判断 出 木 芙蓉 和 台湾 芙蓉 的 杀 缘 
KRR, H UBER 5 CAESUS 最 为 近似 。 一 些 陆 生 植物 ， 如 花生 、 樱 桃 和 葛 卜 等 ， 叶 绿 体 基 因 组 
的 大 小 和 结构 , 以 及 基因 的 含量 和 顺序 , 在 栽培 种 和 野生 近 缘 物种 中 高 度 保守 (Cho et al., 2018; Wang et al., 
2019)。 木 芙蓉 3 个 品种 与 其 近 缘 种 台湾 基 蕉 进化 特征 高 度 保守 ， 再 次 证 实 了 这 一 结论 。 木 芙蓉 与 台湾 芙蓉 
多 态 上 的 差异 非常 小 ， 区 别 仅 为 植株 的 毛 被 类 型 ， 而 木 芙 敬 与 木村 和 朱 检 在 形态 上 则 有 较 大 的 区 别 。 木 英 

车 与 台湾 芙蓉 杂交 坐果 率 较 高 ， 而 与 朱 检 和 木 树 杂 交 坐 果 率 很 低 甚 至 无 法 杂交 成 功 。 张 露 等 (2021) 在 木 
芙 芍 品 种 花粉 的 研究 中 根据 形态 将 “ 单 瓣 白 ” 和 “金秋 颁 ” 聚 在 一 起 ， 而 “牡丹 粉 ” 聚 在 另 一 棱 。 因 此 ， 
本 文 的 研究 结果 验证 了 形态 和 育种 上 的 事实 ， 能 够 较 好 地 研究 进化 特征 和 判断 杀 缘 关系 的 远近 。 此 外 。 木 
芙蓉 3 个 品种 是 一 个 杂交 组 合 ， 其 中 “ 单 激 白 ” 是 “金秋 颁 ” 的 母 本 ， 二 者 最 为 近似 ， 因 此 可 以 判断 出 木 
pe 芙蓉 的 cp DNA 为 母系 遗传 ， 这 为 木 芙蓉 育种 和 遗传 的 研究 打下 了 一 定 的 理论 基础 。 
I~ 3.2 木 芮 蓉 及 近 缘 种 的 叶绿体 基因 组 系统 发 育 关 系 
C J 520 B MUR, ANTEA I — P SRESSAETÉ EIU ERRARE, KRENA Ë 2889272588, 
RARAWA ZBE. KREM F 3 个 品种 聚 成 一 个 支持 率 999085] 8. 28 S0, PA LS 6325 35 2 zÉ pk. 
一 个 高 支持 率 (10090) 43, AEBIAKARAURIG BOSE AER A nut; MRANA, RRR BUM. X 
BRE EXER. AR Ax. BIRRE, GARRIAK: TRÉCOTCNBERUKRTES KER ES SS DN 
Bh RE. AER. DOBRERTEGRARCOKOR EGET. WR Y cp DNA 的 结构 特征 可 以 反映 出 木 芙 敬 3 个 品种 与 
台湾 芙蓉 的 亲缘 关系 外 ， 系 统 发 育 分 析 能 更 直观 的 反映 出 这 一 关系 。 在 系统 发 育 树 中 ， 木 芙 芍 3 个 品种 聚 
在 一 起 ， 与 其 亲缘 关系 最 近 的 为 台湾 美 著 ， 与 木 棋 和 朱 栓 的 杀 缘 关系 则 相对 较 远 。 此 外 ， 建 树 结果 显示 木 
检 属 形成 了 一 个 单 系 。 在 冯 国 棉 (1984) 的 分 类 处 理 中 ， 他 按照 宏观 形态 将 棉 属 和 桐 棉 属 归 入 木 检 族 。 但 
HÆ (2018) 在 息 粉 学 研究 表明 ， 棉 属 、 桐 棉 属 与 木 样 属 的 花粉 有 较 大 差异 ， 与 锅 萎 族 的 一 些 物种 花 
粉 差 异 很 小 。 本 研究 中 , 棉 属 、 桐 棉 属 与 锦 蓉 族 聚 在 一 起 , T] ES IRR AOT» 文 持 之 前 孢 粉 学 的 研究 结果 。 
系统 发 育 分 析 结 果 表 明 叶 绿 体 基因 组 数据 对 于 研究 系统 发 育 问题 是 非常 有 效 的 。 
3.3 木 芮 蓉 及 近 缘 种 的 叶绿体 基因 组 分 子 条 形 码 

从 重复 序列 和 核 昔 酸 多 样 性 分 析 中 可 以 开发 出 候选 分 子 标记 和 DNA 条 形 码 ， 可 以 作为 品种 鉴定 的 分 
子 和 条码。 我们 筛选 出 了 木 栓 属 下 种 级 、 木 芮 敬 种 下 3 个 品种 级 cp DNA 的 候选 特异 性 分 子 标记 和 DNA 条 
攻 码 ， 利 用 所 选 算 择 的 这 些 特异 性 的 标记 ， 可 将 种 级 水 平 物种 区 分 开 ;， 再 将 其 用 于 观赏 花 开 价值 极 高 的 森 
木 美 获 种 下 品种 间 ， 亦 可 将 “ 单 兴 白 ”、 “金秋 颁 ”、 “牡丹 粉 ” 这 3 个 区 分 开 。 

微 卫 星 和 重复 序列 广泛 存在 于 cp DNA 中 , 可 用 于 种 群 遗传 学 研究 和 分 子 标记 的 开发 (Bull et al., 1999; 
Chen et al., 2015)。 木 芙蓉 的 微 卫 星 结构 与 木 梓 和 朱 检 有 明显 的 区 别 ， 但 与 台湾 芙蓉 区 别 不 太 明 显 ， 且 在 木 
英 共 品种 内 部 ， 通 过 微 卫星 结构 无 法 区 分 “ 单 瓣 和 白 ” 和 “金秋 颁 ” 。 相 比较 而 言 ， 重 复 序列 结构 分 析 则 能 
明显 地 区 别 台湾 芙蓉 和 “ 木 芙 蓉 ”， 且 能 对 木 芙 芍 3 个 品种 进行 明显 区 分 。 所 以 ， 重 复 序列 分 析 能 够 更 加 
有 效 地 鉴定 木 美 其 品种 。 叶 绿 体 DNA 条 形 码 是 指 具 有 足够 变异 的 短 DNA 序列 ， 通 常 被 认为 有 助 于 快速 、 
准确 地 鉴定 物种 ( 程 芳 婷 等 ，2015$)。 通 过 核 苷 酸 多 态 性 分 析 我 们 筛选 出 了 台湾 芙 芍 与 木 芙 芍 之 间 的 叶绿体 
DNA 条 形 码 ， 包 括 trnR-UCU-atpA. psbZ-trnG-GCC. infA-rps8 和 ndhE~ndhG。 也 筛选 出 了 木 英 次 3 个 


基因 间隔 区 ， 这 是 因为 基因 间隔 区 处 于 较 弱 的 选择 压力 下 ， 并 且 具 有 比 基 因 更 高 的 进化 速度 ， 它 们 更 适合 
于 低 分 类 阶 元 的 系统 发 育 和 进化 研究 ( 陈 士 林 等 ，2009)。 人 台湾 芙蓉 与 “ 单 淤 白 ” 差 异 最 大 ， 有 7 个 基因 和 
14 个 基因 间隔 区 出 现 不 同 程度 的 差异 。 “牡丹 粉 ” 与 “ 单 办 白 ” 之 间 的 差异 仅 有 一 个 ndhB 基因 和 3 个 基 
因 间 隔 区 。ndhB 基因 属于 NADH 脱 氢 酶 基因 ， 其 缺失 突变 会 导致 光合 作用 中 碳 同化 能 力 严 重 下 降 ， 植 物 


可 以 通过 对 该 类 基因 进行 RNA 编辑 来 保证 光 调 节 作 用 ( 刘 佳 等 , 2019)。 有 研究 证 明 大 豆 的 ndhB 基因 在 水 
稻 中 超 表 达 ， 提 高 了 水 稻 的 光合 效率 ， 进 而 影响 了 水 稻 的 农艺 性 状 〈 王 晓 曼 ，2016)。 因 此 ， 我 们 推断 “ 牡 
丹 粉 ” 花 径 很 大 ， 花 量 更 多 ， 植 株 也 更 高 大 ， 或 许 与 ndhB 基因 的 改变 有 关 。 “金秋 颁 ” 与 “ 单 兴 白 ” 差 
异 最 小 ， 仅 有 一 个 重 白 编码 基因 yfl 基因 发 生 蔡 换 。ycjl 基因 是 具有 未 知 功能 的 开放 阅读 框 ， 但 是 烟草 基 
因 敲 除 实验 证 明 ycfl 基因 编码 对 细胞 存活 至 关 重 要 的 产物 (Drescher et al., 2000)。 与 其 他 叶绿体 基因 相 比 ， 

ycfl 基因 进化 速度 快 , 因此 也 是 植物 中 非常 有 用 的 分 子 标记 , 能 够 揭示 低 分 类 水 平 上 的 系统 发 育 关 系 (Amar， 

2020)。 因 此 , 木 芙蓉 品种 之 间 以 及 与 近 缘 种 之 间 也 可 以 通过 一 些 基 因 间 隔 区 和 个 别 特殊 基因 (如 ycfl. ndhB) 
进行 快速 、 准 确 地 鉴定 。 
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